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AHHOTaU,HH 

A review of the photon-number tomography and symplectic tomography as examples of 
star-product quantization is presented. The classical statistical mechanics is considered within 
the framework of the tomographic representation. 

1 INTRODUCTION 

In quantum mechanics, the state is described by the wave function or density matrix (density 
operator). In classical mechanics, the state of system with fluctuations is described by the probability- 
distribution function. The quantum and classical natures are described using different objects 
- operators and functions. But for more easy understanding the nature of systems which has 
classical and quantum parts, it is necessary to have the same language for both flelds. So, we need 
the representation in which classical and quantum states are described by the same object - the 
probability distribution function. For this purpose, the tomographic probability representation 
of quantum mechanics for continious variables called symplectic tomography was introduced 
in [1-3]. The probability representation of quantum mechanics was extended to the case of discrete 
variables in infinite domain (photon-number tomography) in [3-5]. The name photon- number 
tomography was given in [3] for elucidating the physical meaning of measuring the quantum state 
by means of measuring the number of photons (i.e., photon statistics). Photon-number tomography 
is the method to reconstruct the density operator of quantum state using measurable probability- 
distribution function (photon statistics) called tomogram. The tomographic-probability representation 
of classical states was introduced in [6-9]. In the probability representation, the quantum states 
and the classical states are described by the same objects - tomograms. Tomograms are positive 
measurable probability-distribution functions of random variables. In order to describe observables 
by functions instead of operators in quantum mechanics, the quantization based on star-product 
of functions is used. In [10, 11] it was shown that the symplectic-tomography scheme is a new 
example of quantization based on star-product of functions - symbols of operators. In [12, 13] it 
was shown that photon-number tomography is an example of the star-product quantization. 

The aim of this paper is to present a review of the tomographic approaches and their connection 
with star-product quantization scheme. 
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2 GENERAL STAR-PRODUCT SCHEME 

Following [10,11], we consider a density operator p acting in a given Hilbert space. Let us suppose 
that we have a set of operators W(x) acting in the Hilbert space H ^ where the n-dimensional vector 
X = (xi, X2 . . . , Xn) labels the particular operator in the set. We construct the c- number function 
lyp(x) using the definition 

«;^(x) = Tr(pW(x)). (1) 

The function ifp(x) is called tomogram, it is the symbol of density operator p, and the operators 
U{yi) are dequantizers [14]. We suppose that this relation has the inverse, i.e., there exists the set 
of operators '£'(x) acting in the Hilbert space such that 

p = j u;p(x)P(x) dx. (2) 

The operators li'(x) are quantizers [14]. Formulas ([T]) and ^ are selfconsistent, if one has the 
following property of the quantizers and dequantizers: 



Tr 



W(x):D(x') =(5(x-x'). (3) 



Relations ^ and ([2|) determine the invertable map of the density operator p onto function - 
tomogram Wp(x). We introduce the associative product (star-product) of two tomograms Wp^(x) 
and Wp^iiL) corresponding to two density operators pi and p2, respectively, by the relations 

WpiP2(x) = Wp,(x) * Wf,^{^) = Tr (pip2W(x)). (4) 
The map provides the nonlocal product of two tomograms (star-product) 

Wpi (x) * (x) = J Wp^ {x")wp2 (x')ir(x", x', x) dx.' dx". 

The kernel of the star-product is linear with respect to the dequantizer and nonlinear in the 
quantizer operator 

K(x",x',x) = Tr [P(x"):P(x')W(x)" . 

The associativity condition for tomograms means that the kernel of star-product of tomograms 
ii'(x",x',x) satisfies the nonlinear equation [14] 

J ir(xi,X2,y)ir(y,X3,X4)rfy = y" A"(xi, y, X4)ir(x2, xg, y)rfy. (5) 

3 SYMPLECTIC TOMOGRAPHY 

Now we consider the symplectic-tomography scheme [1] as an example of the star-product quantization, 
following [10,11]. In this scheme, the tomographic symbol of the density operator p called symplectic 
tomogram Wp(x) is obtained due to fomula ([T]) by means of the dequantizer 



W(X, p, v) = S{Xi — pq — fp), 



where vector x = {X, /i, u) has the coordinates which are real numbers and 1 is identity operator. 
The operators q and p are the position and momentum operators, respectively. The symplectic 
tomogram w {X, fi, u) depends on the two extra real parameters n and u. The physical meaning 
of the parameters /i and u is that they describe an ensemble of rotated and scaled reference 
frames in which the position X is measured. For = cos (p and u = sin (p, the symplectic 
tomogram coincides with the marginal distribution for the homodyne-output variable used in 
optical tomography [15,16]. 

Symplectic tomogram can be determined as the expectation value of the delta-function calculated 
with the help of the density operator 

w{X, fi,u) = {6{fiq + up-X)). (6) 
The symplectic tomogram is nonnegative function 

w (X, n, v) > 0, 
and it is normalized with respect to the variable X 

j w (X, /i, v) dX = 1. 

The tomogram is a homogenious function 

w (AX, Xfi, Xv) = \\\^'^w (X, /i, u) . 

The density operator of the state can be reconstructed in the symplectic tomography scheme, in 
view of ([2]), with the help of quantizer 

P(X, fi,u) = — exp (^iXi — iup — i^iq^ . 

Since the density operator determines completely the quantum state of a system and, on the 
other hand, the density operator itself is completely determined by symplectic tomogram, one 
can use for describing quantum states symplectic tomograms (positive and normalized) which are 
probability-distribution functions analogous to the classical ones. The quantum state is given if 
the position probability distribution w (X, /i, u) in an ensemble of rotated and squeezed reference 
frames in phase space is given. The information contained in symplectic tomogram w (X, /x, ly) 
is overcomplete. To determine the quantum state, it is enough to know the values of symplectic 
tomogram of the state for arguments satisfying the following condition: (/i^ + z/^ = 1) , where 
H = cos (p. 

The kernel of star-product of two symplectic tomograms of density operators pi and p2 has 
the following form [10]: 

'^(^(^l + /^2)-K/^l + //2)) 
A(Xi,/Zi,Z/i,X2,/i2,I^2,X^, I^) = ^ — ^ 

X exp (^{ {.,f.2 - .2/.0 + 2X, + 2X2 - ^A!^l±^}) . 



4 PHOTON-NUMBER TOMOGRAPHY 

The photon-number tomogram defined by the relation 

uj{n,a) = {n \ D{a)pD^^{a) \ n) 



(7) 



is the function of integer photon number n and complex number a = Re a + i Im a, where p is the 
state density operator and D{oi) is the Weyl displacement operator 

D{a) — exp(Q;a^ — a* a). 

The photon-number tomogram is a symbol of the density operator 

uj{n, a) = Tr [pW(x)] , (8) 

where dequantizer W(x) reads 

W(x) = b{a)\n){n\b-^{a), x = (n,a). 

The density operator can be reconstructed from the photon-number tomogram with the help of 
the inverse formula [3-5] 



E 

n=0 



A (fa /s _ i\ («++"*)(«+")-^ 



7r(l-s2) Vs + 1 



u){n, a), 



(9) 



where s is an arbitrary ordering parameter [17]. The quantizer in the photon- number-tomography 
scheme is 



P(x) 



TT 1 - S2 



s-l 

S + 1 



{a) +a*){a+a)- 



(10) 



It is known [17] that the Wigner function [18], which corresponds to the density operator p, is 
given by the expression 

W^{q,p) = 2Tr {-if'b^-^)] , 



(11) 



where D(I3) is the Weyl displacement operator with complex argument 

with a and being the photon annihilation and creation operators. 
Let us introduce the displaced density operator 

The Wigner function, which corresponds to the displaced density operator, is of the form 

W^Aq,p) = 2Tr \p^b{l3) {-lf'b{-p) . 



(12) 



(13) 



By inserting the expression for the displaced density operator into (fT3l) . one arrives at 



, rt.t a, 



D 



-1, 



WpM.v) = 2Tr [b-\a)pb{a)b{fi) (- 
In view of the properties of the Weyl displacement operator 

b{P)b{a) = b{l3 + a) exp [i Im {13a*)] , 

b-\a)b-\[3) = (b{p)b{a)y\ 

formula (fT4l) can be simplified as follows: 



(14) 



a] 



D(-a) 



a 



(15) 



One can see that the Wigner function ( 1131) corresponding to the displaced density operator is 
equal to the Wigner function (|2T]1 corresponding to the initial density operator but with displaced 
arguments. The photon-number tomogram is the photon distribution function (the probability to 
have n photons) in the state described by the displaced density operator pa (fT2l) . i.e., 



t^(n,a) = P„(a) = {n\po\n), 



n 



0,1,2,... 



(16) 



As an example, we consider, following [19,20], the photon-number tomogram of the Gaussian state 
of one-mode light described by the Wigner function of generic Gaussian form 



W{q,p) 



detcr(t) 



exp (--Qo--^(t)Q^ 



(17) 



where Q = (^p — (p),'? — and the matrix a{t) is a real symmetric quadrature variance 
matrix. The photon distribution function for one-mode mixed light was obtained explicitly in 
terms of the Hermite polynomials of two variables in [21]. The Hermite polynomials of two 
variables i?^^^(|/i, 2/2), where ni,n2 are nonnegative integers and R is a symmetric 2x2 matrix, 
are determined by the generating function 



exp 



--{xi X2) 

n^ no 
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Applying the scheme of calculations similar to the one used in [21] to our photon-number tomogram (flGll . 
we arrive at the photon- number tomogram as a function of the Hermite polynomial of two variables 



uj{n, a) 



yi{a),y2{a) 



(18) 



where the matrix R, which determines the Hermite polynomial, reads 



R 



1 + 2T + 4rf 



2 (a 



pp 



(Tqq liCTpq) 



Ad 



2 (a.. 



pp 



Oqq + 2iapq) 



pq 



Here d is the determinant of real symmetric quadrature variance matrix cr(t), i.e., d = cfppcrqq — a. 
and T is its trace T = app + agq. The arguments of the Hermite polynomial are 

+ {(Tpp - agq + 2iapq) ((g) + i{p) + V2a)] . (19) 



?/i(") = 2/2 (") 



For the state with displaced Wigner function (fT3l) . the probability to have no photons Poi^a) reads 



Po(a) 



exp 



X exp 

where L = l + 2T + Ad. 



4a 



pq 



L 



{2aqq + 1) ((p) + V2Ima)' + (2(Tpp + 1) ((g) + V2Rea 
((p) + y2Ima) ((g) + 72 Re 



(20) 



5 CLASSICAL STATES 

Let us consider the classical state, for example, the state of a classical particle with one degree of 
freedom with unit mass. We suppose that the particle's position g and momentum p fluctuate due 
to the interaction between the particle and some environment. The particle's state is described 
by a probability distribution function /ci(g,p), which is nonnegative /ci(g,p) > and normalized. 
Here we assumed the normalization condition for /ci(g,p) 

^ J fci{q,p)dqdp = 1, (21) 

analogously to the normalization condition of the Wigner function W{q,p). 

The tomogram of this state is determined by the Radon transform of the probability distribution 
function fci{q,p) 

Wf (X, fi, iy) = J /pj (g, p) 5 {fiq + up - X) dq dp . (22) 

It is called classical tomogram. The classical tomogram can be determined as the expectation 
value of a delta- function calculated with the help of the distribution function /ci(g,p) in the phase 
space 

Wf{X,n,u) = {5{fiq + iyp-X)). (23) 
The classical tomogram is nonnegative function and it is normalized 

J Wf {X, u) dX = 1. 

The physical meaning of the classical tomogram is that it is the probability density for the particle's 
coordinate X, which is measured in the phase-space reference frame subjected to the scaling of 
axes and subsequent rotation with respect to the original reference frame. In the same way as in 
the case of Fourier transformation, where information contained in the function is equivalent to 
information contained in its Fourier transform, information on the particle's state contained in 



the distribution function fd{Q,p) is equivalent to information contained in its Radon transform - 
the classical tomogram Wf {X, /i, v). The Radon transformation is invertible 

/cl P) = ^ / ^/ ^) exp [-1 inq + yp- X)] dX d/i dv, (24) 

so the probability distribution function /ci(^,p) can be reconstructed in view of the classical 
tomogram Wf (X, //, v). 

The commutative star-product kernel K{X, n, u, Xi, ni, i/i, X2, 112-1 ^2) for two classical tomograms 
Wfi{Xx,iii,vi) and Wf2{X2, fJ,2,J^2) is [9] 

Koi^sUX, fx, u, X,,iii, u,,X2, 112, V2) = ■^e^(^i+^^-^('^^+^^)/'')(5(i.(//i + fi2) - fii'^i + 1^2)) . 

The relationship between tomographic star-product kernels in quantum and classical mechanics 
reads [8] 

Equant A*, l^: ^1, t^l, ^1, ^2, ^l2, ^2) = i^classic(^, A*, l^: ^1, t^l, ^2, t^2: U2)e^'^''''''-^''''^^^''l 

6 CONCLUSIONS 

We reviewed two tomographic-probability approaches — symplectic tomography and photon- 
number tomography and showed that they are examples of the star-product quantization scheme. 
We mentioned that, in the quantum case, the state of the system can be described by the 
probability distribution function instead of the density operator. This provides the possibility to 
describe the objects of both the quantum and classical natures using the same objects - tomograms. 
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